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Grover's quantum search algorithm drives a quantum computer from a prepared initial state to a 
desired final state by using selective transformations of these states. Here, we analyze a framework 
when one of the selective trasformations is replaced by a more general unitary transformation. Our 
framework encapsulates several previous generalizations of the Grover's algorithm. We show that 
the general quantum search algorithm can be improved by controlling the transformations through 
an ancilla qubit. As a special case of this improvement, we get a faster quantum algorithm for the 
two-dimensional spatial search. 
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I. INTRODUCTION 

The search problem is to find a desired item sat- 
isfying certain properties out of a given database of 
N items. Consider a quantum system with an A- 
dimensional Hilbert space Hn, whose basis states \j), 
j € {0, 1, ... , A— 1}, encode the A items of the database. 
The target state \t) corresponds to the desired item, while 
all other basis states are non-target states. In general, the 
system can be in any normalized superposition of the ba- 
sis states. We initialize the system in the source state |s) 
and assume a = \(t\s)\ to be non-zero. Direct measure- 
ment of |s) yields the target state \t) with a probability 
a 2 . So on average, 1/a 2 preparations of |s) and subse- 
quent measurements will yield \t) with high probability. 

Grover's algorithm, or more generally quantum am- 
plitude amplification [J, 0, Hj], makes the search faster, 
by performing measurement only after transforming \s) 
to \t) using 7r/4a iterations of the unitary operator Q = 
—I s It, where 1^ = Ijy — 2\tp) is the selective phase in- 
version of the state That is 0(l/|(i|s)|) times faster 
than the measurement scheme. Generally, \s) is prepared 
by applying a unitary operator U on a particular basis 
state, say |0). Then \s) = U\0), and \t) is obtained by 
iterating the operator Q = —UlolI'It on U\0). Grover's 
algorithm has been proved to be optimal, i.e. no other 
algorithm can get to \t) faster than it does 

To search a database using Grover's algorithm, we 
choose \s) to be the equal superposition of all the basis 
states, i.e. |s) = J2j Such a state is generated 

by applying the Walsh-Hadamard transformation W on 
the basis state |0), i.e. |s) = W\0). Then Grover's algo- 
rithm iterates the operator Q = —I s h = —WloWIt on 
\s) to reach \t). If \t) is a unique basis state, a = 
and |t) is reached using [tt^/N /4J queries. 

Grover's algorithm has been generalized in several 
ways. One generalization is to replace Q by Gip,4> = 
e-^Iflf, where 1$ = 1 N - (1 - e lu} ) is the selec- 
tive phase rotation of by angle ui (clearly Tl = 1^). 
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Then \t) is reached by applying 0(csc-^/a) iterations 
of Q<p,$ to |s), provided the phase-matching condition 
tp — 4> — 0(a) is satisfied 0, |(|. Kato's algorithm Q 
replaces I s by an operator K s that, unlike I s , consists 
of only single qubit operators and hence is easier to im- 
plement physically. Ambainis' algorithm 8] replaces I s 
by a real operator i? s , having |s) as its unique eigenstate 
with eigenvalue 1. It has important applications in fast 
quantum walk al gori thms for element distinctness [8J] and 
spatial search [9L llOL [HI ] . 

Here we study further generalizations of Grover's al- 
gorithm. Section II presents and analyses a general 
quantum search algorithm that iterates the search op- 
erator S = D s lf on \s) to reach \t), where D s can be 
almost any unitary operator with |s) as its eigenstate. 
We find the conditions that D s must satisfy to yield a 
successful quantum search algorithm, i.e. reach \t) from 
\s) using 0(1/ a) iterations of S. Section III illustrates 
the general algorithm by explaining Grover's algorithm, 
phase-matching condition, Kato's algorithm, and spatial 
search as its special cases. Section IV presents a con- 
trolled algorithm where the {D s ,lf} operators are con- 
trolled through an ancilla qubit to speed up the search. 
It provides an important application for the two dimen- 
sional spatial search problem, solving it in 0(V N In A) 
time steps [l2j. Earlier algorithms solved this problem 
in 0(y/~N In N) time steps, and it was an open ques- 
tion to design a faster algorithm. Another application is 
to find an eigenstate of a given unitary operator, corre- 
sponding to a known eigenvalue. These results are easily 
extended, by time reversal symmetry, to problems with 
interchanged roles s <-> t. 



II. GENERAL QUANTUM SEARCH 
ALGORITHM 

A. General search operator 

In this section, we analyse a general quantum search 
algorithm that iterates the search operator S = D s lf on 
\s) to take it close to \t). Here D s is a unitary operator 
with the initial state |s) as its eigenstate. We want to 
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derive conditions on the eigenspectrum of D s that allow 
an efficient construction of the search algorithm. Since 
a global phase is irrelevant in quantum dynamics, we 
choose the convention D s \s) = \s). In case D s has an 
M-dimensional degenerate eigenspace with eigenvalue 1 , 
orthonormally spanned by \s m ) with m G {1, 2, . . . , M}, 
we choose the initial state Is) to be 



a 2 = Yj{s m \t)f 



(1) 



Note that a = |(i|s)|. We consider a < 1 so that the 
quantum search algorithm provides a significant speedup, 
0(\/a), over classical search algorithms. 

Let the normalized eigenstates \£) of D s satisfy 

D s \£) =e te '\£), D\\l) =e~ ie *\l), 



0, \e^ s 



> 



>0, 9i G [-tt.tt]. (2) 



the definition Jf 



/flA* 



Here the lower bound m j n characterizes the gap in the 
eigenspectrum of D s . 

To analyse the iteration of S on \s), we need to find 
its eigenvalues and eigenvectors. Let |A/-) be the nor- 
malized eigenvectors of S with eigenvalues e lXk . Now, 
Ijv - (1 - e^)|t)(i| implies that 
|A fc ) - (1 - e^){t\\ k }\t). Left multiplication 
by (£\D S = e l6 '-{£\ gives 

(£\D 8 lf\X k ) = e^{£\X k ) 

= e"' [(£\\ k )-(l-e^)(t\\ k )(£\t)} . 

Thus 

Using (t|A fc ) = EAWiW, we get 



e 

which leads to 



(3) 



(4) 



o i(\ k -0 e ) 



%4> 



^|(l|t)| 2 cot. 



— = cot — 

2 2 



(5) 



The last result follows from (1 — e M ) _1 = | + § cot | 
and the normalization condition |(^|£)| 2 = 1- It is the 
essential condition for e lXk to be an eigenvalue of S. 

Since the R.H.S. of ([5|) is a constant, and cot a; varies 
monotonically with x except for the jump from — oo to 
+oo when x crosses zero, there is a unique solution Afc be- 
tween each pair of consecutive Ogs. Then with 0i =Sm = 0, 
there can be at most two solutions Afc in the interval 
[— 0min, 0min]- It is also clear that one can focus on the 
solutions Afc on either side of any particular Og, by a suit- 
able choice of the inconsequential global phase, i.e. by 
considering the search operator e~ l9e D s lf . 



B. Eigenvalues of the general search operator 

To explicitly solve the secular equation {5| , we assume 
that the smallest eigenphases of S obey 



Assumption : \X k \ -C 6 n 



I A fe | « 1. (6) 



It makes the £ — s m contribution the most significant 
one in the sum on L.H.S. of |JSJ). Now we use the Laurent 
series expansions, ignoring 0(A 2 ) terms, cot 4^ = x~ — if 
and 



Afc — ®t @i X k ( ■> 6t 

cot — - — = — cot — — 1 + cot — 

2 2 2 V 2 



Along with 9 e=Sm = and J2 m \( s ™\t)\ 2 = a 2 , © then 
reduces to a quadratic equation 



2a 2 



d2 

A Afc =0. 

2 



(7) 



Its two solutions are the expected eigenphases on either 
side of 0i= Sm = and satisfying ©. The coefficients 
A, B are, up to 0(a 2 , X 2 



kJ- 



A = cot ■ 



Ax , B 2 = 1 + A 2 , 



where 



A P = Y, I(^)I 2 cof| , pe{l,2}, 



(8) 



(9) 



is the p moment of cot with respect to the distri- 



bute 



over all £ ^ s m . While Ai can have 



any sign, A 2 is always positive. Since 9i G [— 7r,7r], 
we have | cot ^ \ < t^t < j 2 — for £ ^ s m . Then 
E^ Sm IW)| 2 = 1 - a 2 < 1 gives the bounds 



|A P | < (2/e min r, B 2 G [1, 1 + (4/Ci„)] 



(10) 



Moreover, putting x t = \(£\t}\ and y e = \(t\t)\ cot ^ 
in the Cauchy-Schwartz inequality, (^xiyi) 2 < 
(E^)(E2/I) I gives 



A 2 < (1 - a 2 )A 2 < A 2 



(11) 



Let the two roots of ([7]) be A±. Their product is 
A+A_ = —Aa 2 /B 2 so we can write 



2a , . +1 
A± = ± — (tanr?)* 1 , r? G 



0, 



(12) 



The sum of the two roots determines the angle if. We 
have A + + A_ = -2A/B 2 = (2a/B)(ta,mj - cot 77) = 
-(4a/5)cot2?7. Thus 



cot 2i] 



A 
2aB 



(13) 
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For all 77, we have | tan 77 — cot 77I < max± (tan 7y) ±:L < 1 + 
I tan?7 — cot r]\. Together with | tan 77 — cot?i| = |2 cot 2?7| , 
we then obtain 

2\A\/B 2 < max ± |A±| < (2a/ B) + (2\A\/B 2 ). (14) 

Thus the assumption ([6]) for the validity of our analysis 
necessarily holds, when 

2aB + 2\A\ « 9 min B 2 e [0 mln , 9 min + 40jj , (15) 



where \w) is the normalized projection of \t) on the |A±)- 
subspace, and |Aj_) is an unnormalized component or- 
thogonal to this subspace. 

Note that ||A_l|| 2 = 1 - S" 2 csc 2 | > requires 
S- 2 csc 2 | < 1. With B 2 = 1 + A 2 and esc 2 x = 

1 + cot 2 x, it is equivalent to having cot 2 f < A2. That 
trivially holds for <fi — n. In general, <f> and A2 are 
mutually independent, arising from if and D s respec- 
tively, and the constraint has to follow from the assump- 
tion A± <C #min- In our analysis, both |^4| and |Ai| are 



making use of (H)]). With aB > and 9 min < n, we bounded by constant multiples of 9~l n . Therefore 
certainly need \A\ < ^—^ 



20 „ 



In general, the validity condition (|15| involves param- 
eters of the database, D s and if. We will look at several 
special cases in Sections III and IV. 



cot 2 I = A 2 + A(A-2A 1 ) 



= A 



O 



< A2 + O 



B 2 \X, 



C. Eigenvectors of the general search operator 



We calculate the eigenstates |A±), again dropping 
0(a 2 ,A±) terms from the analysis throughout. We 
choose |A±) such that (f|A±) are real and positive. Then 
the normalization condition J2e l(^l^±)| 2 — 1 m © gives 



esc 2 § 
<*1A±> 2 



E 



2 2 a± 

2 CSC 



(16) 



In the sum over £, the £ — s m terms contribute as per 



esc 2 = 4a 2 /A| = B 2 (tanr ? ) =F2 . 



(17) 



To find the sum of < / s m terms, we use the 



assumption 



to get the expansion esc 



2 A±- 



esc" 



l + A±cot^). Then 



I (€|t>| 2 csc : 



2 A±- 



= B Z 



1 + 



A = 



(18) 



where the last inequality follows from A 2 < A2 and (| 14|) . 
Thus we indeed have cot 2 £ < A 2 [l + 0(X±/9 min )], con- 
sistent with the 0(A±/6> m i n ) accuracy of our analysis. 

Putting (|19p in |3]), we obtain the desired eigenvectors, 
|A±>=£^|A±>|£},or 



|A±) 



-A±)/2 



By/1 + (tan 77) 



T2 



£(^csc^V</ 2 \t). 
e 

(21) 

To find the initial state |s) in terms of |A±), we calcu- 
late the projections (s|A±). Using |T]) and Qi= Sm = 0, 
together with (fl~2)) and l|2l"j) . we get 



(*|A±> 



<»»(*- A±)/2 



By/l + (tanTy)^ 2 a 



|(s m |t)| 2 csc 



A. 



-A±)/2 



±- 



y/l + (tan 7/) 



±2 



(22) 



We observe that £ ± |(s|A±)| 2 = 1, and so we do not 
have to worry about component of \s) orthogonal to the 
I A±) -subspace. Explicitly, we have 



Here we have used the inequality | cot '*J m \ < j 2 — and 



the expression B 2 = J2i^ s 



2 esc 2 % which follows 



from (Jll) and © with esc 2 = 1 + cot 2 ^. Ignoring 
0(\±/9 nlin ) contribution in accordance with ([6|), (fl"6]) 
then gives 



(;|a±) 2 



(t\X d 



csc 2 f 



B 2 [l + (tan 77)^ 



By/1 + (tan7i)=F2 



(19) 



Thus 



\w) 



B\ sin f 



■ + |A X }, =sm77|A+)+cos7/|A_), (20) 



\s) =e-^/ 2 [e lA +/ 2 cos77|A + )-e lA -/ 2 sin7;|A_)], (23) 
which yields 

= e -^/2[ e «g'A+ C0Sry |A + ) - e^'A- s i n ry|A_)], (24) 

where q 1 = q + | . 

The results of this Section hold up to 0(a 2 7 A±/6> m j n ), 
and to that extent the |A±)-subspace suffices for our anal- 
ysis. We note that the eigenvalues are determined to 
higher accuracy than the eigenvectors, as is common in 
unitary quantum mechanics. 



D. Algorithm's performance 

The success probability of the algorithm, i.e. obtaining 
\t) upon measuring <S 9 |s), is Pt(q) = \(t\S q \s)\ 2 . From 
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O and (El, 

\{t\S*\s)\ 
With AA = A+ - A_ 
\{t\S«\s)\ = 



sin 2i]\e 



iq'\- 



2P|sinf | 



4a/(Psin2?7), we get 



sin 2?7 
PI sin £ 



2</a 
P sin 2?7 



(25) 



(26) 



As the maximum value of sin a; is 1 for x = tt/2, the 
maximum success probability P m , and the number of it- 
erations q m needed to achieve it, are 



\u)=S^\s), P m = \(t\u)f 



sin 2r) 
B 2 sin 2 1 



9m = 



7T 




ttB sin 2?7 


_A+ - A__ 




4a 



(27) 



We evaluated the eigenvalues A± of S to 0(a 2 ,X±) 
accuracy, as the actual R.H.S. of ([7]) is 0(a 2 ,X±). 
The eigenvalues of S qia are therefore accurate up to 
0(q m a 2 , 9mA 2 -). Since q m — tt/AX = O^Xg 1 ), the 
0(q m \±) terms have the same size as the 0(X±/9 m i n ) 
corrections to the eigenvectors of S. We will see in what 
follows that useful applications of our algorithm have 
5m = @(l/a:), m contrast to the 0(l/a 2 ) query complex- 
ity of classical search algorithms. That makes 0(q m a 2 ) 
terms O(a), and overall validity of the results (I27p is up 
to O(a,X/0 min ). 

The net query complexity of a general algorithm, which 
achieves 0(1) success probability by 0(1/P m ) times 
preparation and measurement of the state |u), is 



Q 



Pu 



tt B 3 sin 2 \ 
4a sin 2rj 



(28) 



Note that quantum amplitude amplification is a faster 
way to take \u) to \t), by applying Q(l/\/P m ) times 
iterations of — I u I t to it. With I u = S Qm I s S~ q ' n , the 
net query complexity then would be Q(q m /VPm) = 
0(P 2 | sin &\/a). But I u requires implementation of I s , 
and the construction of our general algorithm started 
with the hypothesis that we have only the operator D s 
available and not I s . Without quantum amplitude ampli- 
fication, we are restricted to the query complexity (|28|) . 

A particularly interesting case is A = 0, whence (fT"3|) 
gives r] = j. Then 



Pm = 



1 



B 2 sin 2 



ttB 

4a 



ttB 3 , 6 

sm 2 - 

4a 2 



(29) 



With cost; = sin 77 = -4=. (ppjl and give |(w|5 9 |s)| = 
sin(g'AA/2) = sm(2q'a/B). Thus, A = allows us to 
reach the state \w) by iterating 5 on \s), [irB/4:a\ times. 



That maximizes the overlap of \u) with \t), which is not 
possible for A^O. For <j> = tt, we get the best results, 



Pm = 



Ml 2 = -p ' 9" 



7rP 

4a 



7rP 3 

4a 



(30) 



Comparing with the optimal query complexity of 
Grover's algorithm, 7r/4a, we see that our algorithm is 
close to optimal provided B ^> 1. 

We point out that the difference AA between the rel- 
evant eigenvalues governs the rate of rotation from |s) 
towards \t). The search algorithm accomplishes its task 
by flipping the sign of the relative phase between |A + ) 
and |A_) components of the state, and hence needs 
<Zm = 7r/AA iterations. The query complexity of our al- 
gorithm based on the search operator D s is larger than 
q m , because the target state \t) is reached with proba- 
bility less than one. In Section IV, we will show how to 
enhance this probability to P m = 6(1), by controlling 
the operations of our algorithm with an ancilla qubit. 
Those manipulations improve the query complexity to 
Q = nB/Aa, closer to the optimal result of Grover's al- 
gorithm. 

We also note that any practical discrete oracle should 
distinguish the target state sufficiently well from the non- 
target states. Then = 0(1), and indeed = tt gives 
the best results. On the other hand, our analysis is valid 
for any value of (j). 

Another special case is when D s is a real orthogonal 
operator R s . Then its non-real eigenvalues come in com- 
plex conjugate pairs e lSe + — e~ lSe - , and the correspond- 
ing eigenstates satisfy = \£-)*. The cancellations of 
£± contributions, and cot ^ = for 9t = tt, make Ai = 



and A = cot 
2(cotlf+aP) t) ninl 



^ . Then the validity condition (fTS"]) becomes 
P 2 . Also, gives 



sin 2rj = (1 + cot 2 2ij) 



Consequently, 



P, 



-1/2 



2aB/J4a 2 B 2 + cot 2 § 



(31) 



4a 2 



4a 2 B 2 sin 2 § + cos 2 f ' 



ttP 2 



2 x /4a 2 P 2 + cot 2 



(32) 



With (j) = tt, A — cot f = and {P m , Qm, Q} are given 



by (|30|) . This is the case for Ambainis' generalization of 
Grover's algorithm, which leads to the optimal algorithm 
for the element distinctness problem [8j. 

Lastly, consider the situation when \A\ 3> 2aP. Then 
(fl~3| means |cot2ry| ^> 1 and sin2^ w 1 / 1 cot 2?y| = 
2aB/\A\ < 1. As a consequence, using ([28]) and ([2D]). 

7r|A|P 2 sin 2 | ttB \A\ 1 



8a 2 



4a 2aP \(t\w}\ 2 



» Ta • < 33 > 
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where the last inequality follows from B > 1 and 
|(£|w)| 2 < 1. Thus, for a given B, our algorithm is close 
to optimal only when \A\ ^> 2aB, or when rj is not too 
far off from j. In this situation, which includes the par- 
ticular case A = 0, the validity condition (p~5|) simplifies 
and becomes independent of if, 

\A\ > 2aB =► 2a < min B . (34) 

When an estimate of Ai is available, then it may be possi- 
ble to adjust 4> to satisfy this optimality constraint. Even 
for \A\ > 2aB, as long as |A| < |(i|u;)| 2 = B 2 sin 2 |, 
(|3"3"]l gives Q <§; 7r/8a 2 and our algorithm is significantly 
better than classical search algorithms that have query 
complexity 6(l/a 2 ). 

E. Another general search operator 

Time reversal symmetry allows us to extend our anal- 
ysis to the search operator T = IfD t , where D t has 
\t) as its eigenstate with eigenvalue 1. If D t has an M- 
dimensional degenerate eigenspace with eigenvalue 1, or- 
thonormally spanned by \t m ) with m 6 {1,2, . . . , M}, 
then we consider the target state \t) to be 

i*) = \ £<*mi*>i*m} > « 2 = i^i s )i 2 = £ K*-i a )i 2 ■ 

rn m 

The performance of the algorithm is governed by 
\(t\(I£D t )i\s)\ = \{s\(Dll-v)i\t)\. That can be mapped 
to \{t\(D s If ) q \s)\ by the interchanges 

\s)^\t), <}>^-ip, D s -> Dj, t -+ -Ot . (36) 

With specifying the eigenspectrum of D t , 

A = -(cot|+A 1 ), B 2 = l + A 2 , 

A P = E I0»| 2 cotf| , p 6 {l,2}. (37) 

The interchange keeps a = | {t\s)\ invariant, and our anal- 
ysis holds as long as the condition (fl~5|) is satisfied by 
# m in = mhij^trn \0j\- The sign change A —A can be 
accommodated by r\ — > \ — r\, and the expressions for 
{Pm, Qm, Q} hi (|2T)l and [jzg]) are applicable with the only 
change being </i — > — </?. An interesting feature of the 
evolution using the search operator T is that the |A±)- 
subspace almost completely contains the target state |t) 
but not the initial state \s). 

III. SPECIAL CASES 

The general quantum search algorithm presented in 
Section II encapsulates several previous generalizations 
of Grover's algorithm. In this section, we discuss some of 
them as special cases, providing a different perspective 
on their results. 



A. Grover's algorithm 

In Grover's algorithm, the search operator is Q = 
—I s It- To discuss it as a special case of our general algo- 
rithm, we consider the search operator = e~ ltp Iflf 
with selective phase rotations ip and <f> (Q^^ = Q). Then 
D s = e-^If, e l6e = e-^i 1 -^,*) and 6» min = \ip\. With 
6 S = and 0^ s — —ip, ^ gives 

d> ld sin „ tp , . 

A = cot ~ - cot £ = j— 2 — , B 2 = esc 2 £ . (38) 

2 2 S i n £ S in f 2 V 7 

For a fixed B, the constraint (f3"4"]l demands A ^> 2ai? 
for a nearly optimal algorithm. That requires \(f> — ip\ ^> 
4a sin or \tj> — tp\ <C 1 with a -c 1. This is the phase- 
matching condition When it is satisfied, the valid- 
ity condition (TT5)) of our analysis is also satisfied, because 
q«1 and # m in-B = <pcsc ^ > 2. Note that is a spe- 
cial case of T = IfD t as well, so our analysis also holds 
in case of multiple target states, if we consider \t) to be 
given by (|35]) . 

For (f> = ip, the phase-matching condition is certainly 
satisfied and 77 = tt/4. Then JT2J, U2DJ) and (J23J yield 

A± = ±2asin § , 
(A±|i) = (\±\w)=l/V2, 
(X±\s) = ±e- l<l> / 2 e ±las ' m % /V2 . (39) 

The evolution is thus totally confined to the |A±)- 
subspace. Also, (US]), ((27]) and ((5S} provide 

Kt|^| S )|=sin[2a(g + i)sinM] 5 

Pm = 1, = 9m = k CSC M/ 4a j . (40) 

For tf> = (p = 7r, we have Grover's algorithm. Then 
= sin[(2<? + l)a] and Q = [Tr/Aa\. 

B. Kato's algorithm 

Let the Af-dimensional Hilbert space be spanned by 
n = log 2 N qubits, with each basis state \j) correspond- 
ing to an n-bit binary string j € {0, 1}™. Then the op- 
erator I s cannot be written as a tensor product of sin- 
gle qubit operators. Its implementation needs coupling 
among qubits, which may be difficult to realize physi- 
cally. Kato's algorithm [7] replaces I s by K s , which is 
a tensor product of single qubit operators. The algo- 
rithm iterates the operator K s lf on \s) = H® n \0), with 
K s = (HZjH)® 71 , where 

*=^(i -0' z ^^( e o A)' (41) 

and the phase 76 [— — , ^] as determined below. 
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The eigenvectors \£) of K s are \i) = H® n \j), with the 
corresponding eigenphases 6( — (ir — 2^hj) mo ^ 27! — 7r, 
where /ij € {0, . . . , n} is the Hamming weight of the bi- 
nary string j. For the binary string of all zeroes, ho = 0, 
and so 9i =s — 0. Then the restriction | '-y j < -^-^ en- 
sures m i n = 2|7|. The degeneracy of 0£ is C n ,hji the 
choose function for selecting hj items out of n. Also, 
\{£\t}\ 2 = 2~ n for all L With all these values, © gives 

A p = (-!) P £ 2""^,^ cot p 7 ^ ,pE{l, 2}. (42) 

The binomial weight 2~ n C n j lj is maximum at hj = 
and for large n decays exponentially away from it as 
y/2/im exp[— 2h 2 Jn\, with hj — hj — ^. The signifi- 
cant terms in the sum over hj therefore correspond to 
hj = 0(y/n) <§C ^ for large n. To evaluate the sum in 
(|42|) . we use second order Taylor expansion of cot p jhj 
around 771/ 2, 



(-l)f A p « cot p ^ + 



2 

7 n 



d 2 cot p 7/ 
dy 2 



cot^ 



7 2 n 9 771 
- — esc — 
4 2 



771 



(2p-l)cof -^ + 4 



'p,2 



where we have used the moments ^2 h .2 n C rl ^h j h r ^ = S r 
with So = li Si = and S2 — n/i. Using the bound 
cot |^| < esc |3? I < |^|, we get 



-Ai 



cot — + 0(7) esc 2 — 
2 K " 2 



cot 



jn[l - 0(l/n)] 



A 2 = [1 + O (1/n)] cot 2 ^ + O (1/n) . (43) 



Kato chose 7 such that A = cot $ + Ai = 0. Then ([43 
gives 



7=(<£/n)[l + 0(l/n)], 
A a = [1 + 0(l/n)] cot 2 f + 0(l/n). 



(44) 



Thus B 2 = 1 + A 2 = [1 + 0(l/r7)] esc 2 § + 0(l/n), 
and the validity condition (fT5|) becomes 2a <C # m in-B ~ 
27 esc ^ < 27r/ri for large 77. It is certainly satisfied be- 
cause a = \(s\t) \ — 2~ n / 2 . The performance of the algo- 
rithm follows from (|29|) . 



7T2" 



p m = i-o(i) , 

(1 + 0(1)} cscM + o(i) 



(45) 



Comparing with (|4"D|) , we see that Kato's algorithm is 
almost as efficient as Grover's algorithm for large 77. We 
also have the proof that its maximum success probability 
behaves as 1 — 0(l/n), which was shown only numerically 
by Kato. We point out that the above analysis does not 
hold for <p — tt, due to the restriction I7I < that 

demands \<f>\ = n| 7 |(l - 0(±)) < tt(1 - 6(±)). " 



C. Two-dimensional spatial search 

Consider a two-dimensional square lattice, whose sites 
encode the basis states of an A-dimensional Hilbert space 
TLn- Let the sites be labeled by their coordinates as \x, y) 
for x, y G {0, . . . , VN — 1}, and let \x t ,yt) be the label 
of the target state. A local operator L is an operator 
that does not couple any pair of non-neighboring sites, 
i.e. (x',y'\L\x,y) = if \x' - x\ + \y' - y\ > 1. The 
two-dimensional spatial search problem is to find \xt, yt), 
with the constraints that in one time step we can execute 
either an oracle query or a local operation. Such a situ- 
ation can arise when the database is spread over several 
distinct locations, and locality of physical interactions 
prevents large jumps between locations. 

In Grover's algorithm, the initial state \s) is a uniform 
superposition of all lattice sites, i.e. \s) = ^7= J2 X y \ x > v)- 
To prepare it, we start from a particular site, say |0, 0), 
and then repeatedly transfer amplitudes to neighboring 
sites using a local operator. As we must transfer am- 
plitude across all y/N sites in both directions of the lat- 
tice, the preparation takes 2y/N local operations. Thus 
if \ s ) = W\0,0) then W takes &(Vn) time steps, and 
the same is true for I s = W/o.oWt Since Grover's algo- 
rithm requires one application of I s per query, its total 
time complexity becomes 0{y/~N x \/N) — 0(N), which 
is no better than that of classical algorithms. 

To obtain quantum speedup for spatial search prob- 
lems, Grover's algorithm has been modified in several 
ways [1, |T(J HH • We look at the particular modification 
due to Ambainis, Kempe, and Rivosh (AKR) [9(. Their 
algorithm uses a 4-dim ancilla coin space 7i c , whose basis 
states I — ►), I <—),] 1), I |) represent the four possible di- 
rections of movement on a 2-dim lattice. Let \u c ) be the 
uniform superposition of these coin basis states. AKR's 
algorithm is implemented in the 4A r -dim Hilbert space 
TLc^'Hn- It drives the initial state \s') = \u c )\s) towards 
the effective target state \t') = \u c )\x tl yt) by iterating 
the operator L s il t i. Here I t > = I4JV — 2\t')(t'\ is imple- 
mented using an oracle query, while L a > = CI Ua is a local 
operation with I Uc = 2|u c )(u c | — I4 and C defined by 

\^)®\x,y) -^\^-)®\x + l,y), 
I <-)®\ x ,y) ®\x-l,y}, 
I T)<s>M I l)®\x,y + l), 
I l)®\x,y) I T>®|a;,y-1>. 
AKR found the eigenstates of L s > to be 



(46) 



\£a,b) = \v a fi)\F a )\n), a, b G {0, 1, . . . , VN - 1}, (47) 



where \T a ) = N 4 J2 X e ^ \x) and \F b ) = 
N~i ^2 y e vTv |y) form the Fourier basis. For each (a, 6), 
there are four eigenvalues 1, —1, and e a - b with 



2cos6» Qib = cos(27ra/V7V) + cos(2nb / V N) , (48) 
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corresponding to four orthogonal vectors \v\ 



bh 



and 



u a,b, 



of the coin space. 



AKR have shown that L S 'Iv preserves the subspace 
Ho spanned by the eigenstates |s') = \u c ) |JF ) and 



\F a )\Fb) with (a, b) ^ (0,0). Further- 



more, the states \v a b ) are such that \(£ a b \t')\ — l/\f2N. 
Within Ho, \s'} is a unique eigenstate with eigenvalue 1, 
and L s i is a real orthogonal operator with A = as de- 
scribed in Section II. D. To compute B, we put | (£^ b \t') \ — 

1/V2N and 1 + cot 2 x = esc 2 x = 2/(1 - cos 2x) in (0, 
to get 



B 2 = 



2 1 

TV f-^ 1 - COS 6 a I 

(o,b)54(0,0) 



1 F 

AT ^ 



(o,6)54(0,0) 



27ra 2nb 
2 — cos — ^= — cos —== 

N VN 



(49) 



Using 1 — cosx — 2 sin 2 |, and x > sin a; > ^x for x S 



[0,f] 



E 



(0,6)54(0,0) 



• 2 ^ 

sin" — ^= + sin — == 
N VN 



6 E ?r 

\(q,6)#(0,0) 



/} 2 



e(lniV). (50) 



Here the last equality follows from X)b( fl2 + 
b 2 )' 1 9(l/o)(l « Ol0 ) + e(l)J O|0 , and 

0(Vo) = ©(hiiV). Also, O implies 
that min = min (a:fc )^( 0i0) 9 a<h = ^/2tt/^/N. The 
condition (|15[) for the validity of our analy- 
sis is satisfied for large N, because A = and 
2a = 2/VN < min B = y/QQnN)/N. Thereafter, with 
a = 1/VN and \w') = (L s ,I t ,) q ™\s'), ([3D]) gives 



(i>') 2 = 9(1/ In AT) , ?m = e(VJVlnJV). (51) 



Note that our analysis gives a more accurate estimate of 
P m compared to AKR's estimate, P m = 0(1/ In AT). 

Subsequently, AKR's algorithm uses quantum am- 
plitude amplification algorithm to reach \t') by apply- 
ing 6(1/1 (t'\w') |) = e(VhiA0 iterations of -I w d v on 
\w'}, with I„, = {L s ,I v y™I s .{L\,lly™. Now /,/ can 
be implemented in 0(\/AF) time steps just like I s , so 



I w > takes 6(<7 m ) + 8(vJV) = 6 (VAT hi AT) time steps. 
The net query complexity of the algorithm is, therefore, 
(VAT In AT x VhiA) = G(VFlnAT). 

It was an open question whether or not the 
©(•v/iVlnAT) algorithm for 2-dim spatial search can be 
improved further. In the next Section, we give a posi- 
tive answer to this question by presenting a 0(V N In N) 
algorithm, which is 0(vm7V) times faster than AKR's 
algorithm. 



IV. CONTROLLED QUANTUM SEARCH 

In Section II, we analysed the performance of the gen- 
eral search operator S = D s lf . In this Section, we show 
that just by controlling the operators {D s ,lf} through 
an ancilla qubit, without modifying them at all, we can 
manipulate the moments {Ai,A2} and hence the coeffi- 
cients {A, B). The results $Z7§ and (J25J) imply that the 
best performance of the general quantum search algo- 
rithm is achieved by <fi = it, A = and B > 1 as small 
as possible. The aim of our manipulations is to achieve 
that situation. 



A. Manipulation of the first moment 

To control the search process, we attach an ancilla 
qubit with the Hilbert space Hi to the search space Hn, 
and work in the joint space H' = H2 <8> Hn ■ The logic 
circuit of Algorithm 1 that effectively makes A vanish 
is shown in Fig. 1, where the dashed box represents the 
search operator S' = D' s I t > with 

D' s = ( Cl Dt)(coD s ), 
c D s = \0)(0\®D s + \l)(l\®l N , 
Cl D\ = |0)(0|®1 JV + |1)<1|®£>| ! 
I v = t 2N -2\if){if\, |f>=l+>l*>, (52) 

where |+) = ^|[|0) + |1)]. Here c D s and c\D\ control 

the operation of D s or D\ on Hn through the ancilla 
qubit. The eigenstates of D' s are \£' ) = \0)\£) and ]£[) = 
\1)\£), with 0£ and — 9g as the corresponding eigenphases. 
Thus 

0' to = B t , 9' ei = -9 t , (£' \t') = (£[\t') = -^{£\t). (53) 

D' s has two mutually orthogonal eigenstates, \s' ) = \0)\s) 
and \s[) = |l)|s), with eigenvalue 1. In accordance with 
([I]), we choose the initial state \s') to be 

k') = £ [Wo\«)Wo) + Wi\*)Wi)] = \+)\s)> 

ol = V\(WW+W¥W=\(s\t)\ = a. (54) 

We also have 9' min = min^^ 0# — m in- Now, putting 
(PI in ©, we get 



A;= £ \W)\ 2 «*?f + £ K^i|0| 2 cot^ 



(55) 



Explicitly, A' = cot f + A' : = for 



7T, and B' 



v/TTA^ = B. Since {9' mia , a',B'} = {9 min , a, B}, when- 
ever D s satisfies the condition 9 m [ n B 3> 2a, its equivalent 
condition 9' min B' ^> 2a' is satisfied by D' s . The results 
of Section II hold, therefore, for analysing the iteration 
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l+> — 



Iterate Q' times 

— 1 



(a) 



D s 









— 1+> 
—\t) 



FIG. 1: Logic circuit diagram for Algorithm 1. 



|0> — 
l«> — 



Iterate Q" times 

o 



D s 



— 10 

— 1*> 



(b) 



FIG. 2: Logic circuit diagram for Algorithm 2. 



of D' s I t > on \s') to reach the effective target state \t'), 
whose projection on Hn is the actual target state \t). As 
A' = 0, the corresponding {P^, q' m , Q'} are given by (f2"9")l . 
i.e. 



P' = B~ 



\irB/Aa\ , Q' = ttP 3 / 4 " . (56) 



The importance of this manipulation is that A' is always 
zero irrespective of the value of A, so that we get the 
optimal algorithm for fixed {# m i n , a, B}. 



B. Manipulation of the second moment 

The query complexity Q' can be further improved from 
7rP 3 /4a to irB/4a by manipulating the second moment 
A 2 . Without loss of generality, we analyse only the Ai = 
case, obtained if necessary by addition of an ancilla 
qubit as in Algorithm 1. The logic circuit of Algorithm 
2 that effectively reduces B is shown in Fig. 2, where the 
dashed box represents the search operator S" — D"I t '> 
in the joint space 7i" = H.2 <8 H-n, with 

D% = (Z ® l N )(c D s ), 
\t") = I C> 1*> , 10 =sinC|0)+cosC|l). (57) 

Here Z = |0)(0| - |1)(1| is the Pauli Z operator. The 
eigenstates of D" s are %) = |0)|£) and \£'{) = \l)\i), with 
9t and ir as the corresponding eigenphases. Thus 



(f '|O = (t\t) sinC , {l'{\n = (t\t) cosC ■ (58) 

Here, D" has a unique eigenstate \s'q) = \0)\s) with 
eigenvalue 1, and we choose it as the initial state \s"). So 
a" = (s"|0 = asinC and Cin = min r , #s » |^'| = min . 
Now, putting (J58J) in ©, we get 



K= E i<4'ini 2 cot^ 



IWIO| 2 cot^ 



A" 

P 



A p sin 2 C ■ (59) 

With = 7T and Ai = 0, we have A" = A" = 0. Further- 
more, 



(B") = 1 + A' 2 ' = l + A 2 sin"C 
= cos 2 C + B 2 sin 2 ( C « 1 



For B > 1, we choose C < -S -1 < 1, to obtain 1 < 
(_B") 2 < 2. Also a" w a< < a/S means that 2a/B > 
2a" > 2a" /B", and so whenever £> s satisfies m i n B ^> 
2a, D' s ' satisfies Cin 5 " > 2a "- The analysis of Section 
II therefore holds for D", and (l56l) reduces to 



P" 



1 



(p») 



"1 




TlB" 




, C = 


_4a(_ 



7T (1+P 2 C 2 ) 3/2 

4^ C 



(60) 



Minimizing Q" over C, the optimal choice is C = l/y2B, 
yielding 



i mm 



Q'U 



(3v / 3/2)( 7 rP/4a) . (61) 



1 + P 2 C 2 . 



Thus we have an 0(B/a) search algorithm; the 0(B 2 ) 
speedup is useful only when B 3> 1 that we assumed. 
Note that D\ is not needed here unlike in Algorithm 1. 



C. Applications 

1. Faster algorithm for 2-dim spatial search 

The manipulations of Sections IV. A and IV. B find an 
important application in two-dimensional spatial search, 
discussed in Section III.E. The spatial search is restricted 
by nearest neighbour movement, and in d dimensions 
niN 1 ^) steps are needed to cover the whole database. 
Grover's optimal search algorithm allows movement from 
any site to any other site in just one step, and can be ef- 
fectively considered the d — > oo limit of spatial search. 
These considerations imply that the spatial search prob- 
lem has complexity f2(max(V 1 / d , \/~N)). 

For the d = 1 case, quantum spatial search has the 
same complexity as classical search. For d > 2, one can 
achieve 0(y/~N) complexity with quantum spatial search, 
comparable to Grover's algorithm. The critical d = 2 
case is special, with the conflict between different dy- 
namical featur es pr oducing lnA^ factors. We have seen 
that B = 8(\/ln N) for d = 2, as a result of the diver- 
gence in KSDJ}. (Note that B = 6(1) for d > 2.) With 
a = 1/ V^V, the manipulations described above give an 
O(Bfa) = O(VnToN) algorithm, which is 0(\/In N) 
times faster than AKR's algorithm. 
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2. Finding eigenstate for a given eigenvalue 

Consider the problem of finding an eigenstate \s) of a 
given operator D with known eigenvalue e l6s . We first 
construct the operator D s = e~ l6s D, which has eigen- 
value 1 for \s). The search algorithm then can take us 
to the state \t) using D e , D\ and /(. Reversing the al- 
gorithm, we can start with a known state \t), and evolve 
to \s) using O(Bfa) applications of {I t ,D,D^} (we can 
use Algorithms 1 and 2 if necessary). The required prop- 
erty is that the eigenphases 6e^t s of D should be well- 
separated from 9 S according to the condition Ijl5|> . i.e. 
^min = min^ s (^ - a ) » 2a/ B. 

V. CONCLUSION 

We have analysed a general framework of the quantum 
search algorithms, where one of the selective transforma- 
tions gets replaced by a more general unitary transfor- 
mation. We have derived the conditions for a successul 
quantum search and calculated the number of iterations 
required by the algorithm. We have discussed several 
quantum search algorithms as special cases of our gen- 
eral framework. There are other search algorithms also, 
not discussed here, which can be considered as special 



cases. For example, the quantum random walk search 
algorithm presented by Shenvi, Kempe and Whaley [HI 
is equivalent to spatial search in d = log 2 N dimensions. 

We have shown that the search operators can be con- 
trolled through an ancilla qubit to get faster quantum 
search algorithms. These algorithms may find interest- 
ing applications. For example, a faster quantum walk 
algorithm for the two-dimensional spatial search can be 
obtained using similar techniques [12]. Some other pos- 
sibilities are under investigation. 

We point out that the general framework presented 
here applies only to iterative search algorithms and does 
not apply to recursive quantum search algorithms (see, 
for example, section III of [lij]]). Our analysis provides 
insights in to the nature of iterative quantum search al- 
gorithms, where the performance depends completely on 
the eigenspectrum of the search operator. That is unlike 
the recursive case, where the performance depends upon 
the amplification factor provided by the search opera- 
tor. We believe that our general framework will serve as 
an important tool in designing future iterative quantum 
search algorithms. 
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